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RELIABILITY  ASSESSMENT  FOR  HIGHLY  RELIABLE  SYSTEMS 
by  M.  V.  Johns,  Jr. 

Introduction.  The  purpose  of  this  study  is  to  develop  appropriate 
methods  for  the  assessment  of  reliability  for  highly  reliable  systems, 
i.e.,  systems  characterized  by  failure  probabilities  no  larger  than  5%. 
Particular  emphasis  is  given  to  cases  where  the  probability  of  system 
failure  during  any  particular  duty  cycle  is  less  than  1%  and  where  sub- 
stantial amounts' of  subsystem  test  data  are  available.  Serial  systems 
consisting  of  several  component  subsystems  are  considered  and  confidence 
bounds  for  the  mean  number  of  trials  between  failures  based  on  subsystem 
failure  data  are  developed.  (The  mean  number  of  trials  between  failures 
(MTBF)  is  equal  to  l/(l-R)  when  R is  the  system  reliability,  i.e.,  the 
probability  of  completion  of  a failure-free  duty  cycle. ) Tables  needed 
for  the  implementation  of  the  proposed  procedure  are  provided.  It  is 
assumed  that  the  failure  data  available  for  the  component  subsystems 
consists  of  the  numbers  of  failures,  observed  in  a substantial  number  of 
independent  trials  for  each,  subsystem  where  the  numbers  of  trials  may 
differ  for  different  subsystems.  Previously  available  methods  for  obtaining 
confidence  bounds  for  series  system  reliability  have  typically  required 
the  assumption  of  exponentially  distributed  times  to  failure  (see  e.g.  [3], 

[ 10 ] ) or  have  involved  ad  hoc  approximations  (see  e.g.  [ll]).  Generally 
the  data  are  assumed  to  have  been  obtained  through  inverse  sampling  (Type  II 
censoring),  i.e.,  subsystem  testing  continues  until  a predetermined  number 
of  failures  has  been  observed — a procedure  often  difficult  to  implement  in 
practice  (see  e.g.  [3],  [7J)«  The  approach  used  in  the  present  study  is 
non-Bayesian,  thereby  avoiding  the  difficulties  associated  with  subjective 


probabilities.  The  underlying  rationale  is  based  on  exact  methods,  although 
asymptotic  approximations  are  proposed  and  evaluated. 

Consider  a system  consisting  of  k subsystems  operating  independently 
in  a series  configuration.  Let  p be  the  probability  that  the  ith  sub- 
system will  fail  on  any  given  trial  (duty  cycle)  and  suppose  that  the  sub- 
system is  subjected  to  n trials,  i = l,2,...,k  . We  envisage  a situa- 
tion where  the  pi's  are  substantially  less  than  1%  and  the  n.'s  are  of 
the  order  of  several  hundred  to  several  thousand.  For  such  cases,  if  X 

i 

is  the  number  of  observed  failures  for  the  ith  subsystem,  then  the  repre- 
sentation of  the  distribution  of  X^  by  the  Poisson  approximation  with, 
parameter  X±  = n^  is  essentially  exact.  The  system  reliability  is 
given  by 


k 

R = n tl-p  ) 

i=l 


For  the  values  of  the  p^'s  considered  here  the  approximation 

k 

R = 1 - 7 X. /n. 

i=l  1 1 


is  essentially  exact  so  that  to  find  a lower  confidence  bound  for  the  mean 

time  between  failures  = MTEF  = l/(l-R),  or  for  the  system  reliability,  R , 

k 

it  suffices  to  obtain  an  upper  confidence  bound  for  £ A./n.  . 

i=l  i 

The  problem  thus  reduces  to  that  of  finding  an  upper  confidence  bound 
for  a linear  function  of  Poisson  parameters  based  on  independent  Poisson 
observations  X^  jXg , . . . jX^.  . To  treat  this  problem  we  introduce  the  param- 
eter 


(1) 


k 

0=  l 


i=l 


a.  X . 
i l 


5 
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k 

where  the  a.'s  are  normalized  coefficients  with  Y a.  = 1 and 
1 i=l  x 

ai  > 0 , i = 1,2, ... ,k  . The  use  of  a ’s  which  sum  to  one  facilitates 

the  construction  of  numerical  tables.  The  connection  with  the  original 

problem  is  as  follows:  For  i = l,2,...,k  let 


(2) 


a.  = 

l 


If  0 is  an  upper  confidence  bound  at  confidence  level  1-a  for  0 given 

k 

by  (l),  then  the  corresponding  upper  confidence  bound  for  Y A./n.  is 

» / V \ i=l  1 1 

given  by  0 M l/ru  J , and  the  1-a  level  confidence  statement  for  the 

quantity  MTBF  is 
(3)  MTBF  > 


Because  of  the  additive  property  of  the  Poisson  distribution,  the  data 
for  subsystems  for  which  the  sample  sizes  n..  are  the  same  may  be  combined 
so  that  k represents  the  number  of  different  subsystem  sample  sizes  or 
the  number  of  distinct  a^'s  . If  all  subsystems  are  subjected  to  the  same 
number  of  trials  we  may  take  k equal  to  one  and  the  problem  reduces  to 
the  familiar  one  of  finding  an  upper  confidence  bound  for  a single  Poisson 
parameter . 

The  subsystem  failure  data  may  be  developed  through  independent  testing 
of  the  subsystems,  or  through  testing  of  the  complete  system  with  the  assign- 
ment of  failures  to  the  appropriate  subsystems.  Even  in  the  latter  case 
subsystem  sample  sizes  may  differ  if  subsystems  are  redesigned  during  the 
course  of  testing  so  that  the  trials  and  failures  observed  prior  to  redesign 
for  such  subsystems  are  not  relevant  to  the  reliability  of  the  final  version 
of  the  system. 
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It  is  worth  noting  that  the  3,81116  formal  confidence  "bound  problem  anises 
from  a somewhat  different  and  more  restrictive  model  involving  continuous 
time  sampling  and  exponentially-  distributed  failure  times.  Thus,  if  we 
assume  that  the  time  to  failure  of  the  ith  subsystem  has  the  exponential 

-v 

probability  density  function  y±e  , y > 0 , t > 0 , and  if  the  fixed 

aggregate  testing  time  for  the  ith"  subsystem  is  given  by  t and  the  total 

number  of  failures  observed  is  N.  , i = l,2,...,k  , then  the  N.’s  which 

1 1 

are  sufficient  statistics  have  independent  Poisson  distributions  with 
corresponding  parameters  = y T . Such  test  data  would  be  generated, 
for  example,  if  a single  unit  of  each  subsystem  were  placed  on  test  and 
repaired  whenever  a failure  occurred  until  the  total  test  time  for  each 
subsystem  reached  its  corresponding  limit,  . The  system  reliability 
for  a single  duty  cycle  (consisting  of 'operation  for  a unit  interval  of  time) 
is  given  by 

k -y  / k 

R = n e = exp  - l y 
i=l  \ i=l  1 

If  the  system  is  highly  reliable  the  y^'s  will  be  very  small  and  we  have 


1 - R 


k 

* I 

i=l 


yi 


k 

= l 

i=l 


VTi 


and  the  confidence  bound  problem  becomes  formally  the  same  as  the  one 
previously  introduced. 


- k - 


In  Section  2 the  general  method  for  constructing  the  required  confi- 


dence bounds  is  described  and  optimality  properties  are  discussed.  The 
method  is  based  on  the  use  of  an  appropriate  linear  ordering  imposed  upon 
the  collection  of  possible  sample  vectors  X = (X. ,X_, . . . ,X,  ) and  exploits 
an  idea  first  introduced  by  Buehler  [2]  • The  use  of  the  tables  which  have 
been  computed  for  the  cases  k * 1 and  k - 2 is  illustrated  in  this 
section. 

In  Section  3 approximations  are  introduced  which  permit  the  calcula- 
tion of  the  required  confidence  bounds  for  data  beyond  the  range  of  the 
computed  tables.  These  approximations,  which  are  based  on  the  maximum 
likelihood  ratio  (MLR)  statistic  and  the  asymptotically  equivalent  maximum 
likelihood  estimate  (MLE)  confidence  bound,  are  used  to  justify  a procedure 
for  reducing  cases  where  k •>  3 to  the  k = 2 case.  The  asymptotic 
expressions  obtained  are  not,  strictly  speaking,  "large  sample"  results 
since  it  is  not  the  sample  sizes  (the  n^’s)  which  are  assumed  to  become 
large  but  rather  the  largest  of  the  1 s . 

The  formal  proofs  of  the  optimality  results  of  Section  2 and  the 
asymptotic  results  of  Section  3 are  deferred  to  Appendices  1 and  2. 

Appendix  3 describes  the  construction  of  the  tables. 
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2.  Confidence  Bounds.  A frequent  practice  in  dealing  with  confidence 
bound  problems  is  to  use  procedures,  which  are  "optimal"  in  the  sense  that 
they  correspond  to  a hypothesis  testing  procedure  which  has  some  desirable 
property  such  as  being,  say,  uniformly  most  powerful  unbiased.  Such 
confidence  bounds  are  said  to  be  "uniformly  most  accurate  unbiased"  (see 
[5],  pp.  176-180).  This  approach  is  relevant  only  when  UMP  tests  exist 
and  is  somewhat  questionable  in  any  case,  since  the  criteria  of  interest 
in  confidence  statements  such  as  the  expected  size  of  the  confidence  region 
have  no  counterparts  in  the  hypothesis  testing  context.  For  the  present 
case,  no  UMP  unbiased  or  invariant  tests  of  hypotheses  concerning  linear 
functions  of  Poisson  parameters  exist  so  that  another  approach  is  required. 

As  noted  in  Section  1,  we  have  elected  to  exploit  an  idea  introduced  by 
Beuhler  [2]  which  involves  the  construction  of  confidence  bounds  based  on 
suitably  generated  ordering  of  the  sample  data  vectors. 

Before  discussing  the  specific  problem  of  constructing  an  upper 
confidence  bound  for  the  parameter  0 defined  by  (l),  we  consider  the 
following  more  general  situation:  Let  P^  represent  a family  of  probability 

functions  indexed  by  a (vector)  parameter  X lying  in  a set  A , and 
suppose  that  for  each  X_  , P^  is  defined  on  the  appropriate  subsets  of  a 
sample  space  whose  (possibly  vector)  generic  point  is  designated  by 

x with  X denoting  the  corresponding  random  vector.  Let  0 = &(X) 
represent  a real  valued  function  of  X_  . Suppose  further  that  a linear 
ordering  (denoted  by  ^ ) of  the  elements  of  the  sample  space  is  specified 
and  that  this  ordering  may  be  thought  of  as  being  generated  by  a real 
valued  function  <j>(x)  . I.e.  , if  and  only  if  <j)(x)  >_  <|>(y)  . We 

seek  an  upper  1-a  level  confidence  bound  0*(x)  for  0 which  is  monotone 
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in  the  specified  ordering  on  . Thus  we  require 


U)  PAie  <_  0*(x)>  >_  1 - a for  all  X_  e A , and 

(5)  x <==>  9*W  > e*(y)  • 

To  achieve  this  we  proceed  as  follows:  For  each  real  t in  the  range  of 

0(A)  let  S(t ) = {A;  0(A)  = t>  . 

Assumption  A.  For  each  x e the  function 

h (t)  = sup  P.  {X -<  x} 

- XeS(t)  - 

is  monotone  decreasing  in  t , and  there  exists  a value  t (x)  such  that 

CL 

and  the  collection  of  such,  values  contains  its  greatest  lower  "bound. 

Def.  For  each  x e 3f.  0*(x)  = smallest  t^x)  . 

The  properties  of  0*(x)  are  outlined  in  the  following  four  proposi- 
tions, the  proofs  of  which  are  deferred  to  Appendix  1: 

Proposition  1.  Assumption  A implies  that  0*(x)  is  monotone  in  the 
ordering  (>-)  . 

Proposition  2.  Under  Assumption  A 0*(x)  is  a 1-a  level  confidence 
hound  for  0(A_)  . 

Proposition  3.  Under  Assumption  A,  if  0 (x ) is  any  other  1-a  level 
confidence  hound  for  0(A)  which  is  monotone  in  ( >)  then  0(x)  >_  0*(x) 


Proposition  b.  Under  Assumption  A,  0* (x)  is  undominated  in  the 

sense  that  if  0 (x)  is  any  Itch  confidence  hound,  then  sup  0(y)  >_  0*(x) 

yAx  ~ 

for  all  x_  e % . 

Proposition  3 guarantees  that  0*  is  optimal  with  respect  to  the 
given  ordering  and  Proposition  H insures  weak  admissibility  with  respect 
to  any  ordering,  i.e.,  there  does  not  exist  a 1-a  level  confidence  bound 
0 , say,  such  that  0(x)  < 0*  (x ) for  all  x_  e . Furthermore,  if  for 
every  distinct  x^,  x^  £ > 0*(x^)  ^ 0*(Sg)  » then  for  any  1-a  level 

confidence  bound  0 , if  0 (x ’ ) < 0*(x_’  ) for  some  x1  £ % , then 
0 (x" ) > 0*(x")  for  some  x"  £ ;£  . The  choice  of  the  ordering  (^)  or, 
equivalently,  the  order  generating  function  <j>(x)  remains  to  be  determined. 
In  the  present  application  <£  will  be  chosen  so  as  to  insure  that  0*(x) 
is  reasonable  for  small  true  values  of  0(X_)  and  asymptotically  optimal 
in  the  usual  sense  (uniformly  most  accurate)  for  large  0(A). 


Returning  to  the  specific  problem  under  consideration,  we  recall  that 

k 

X = (X1,X2,...  ,Xk)  , 0(A)  = J a±X±  and  X=  (^  ,Xg , . . . ,Xfc)  where  the 

X.'s  are  independent  observations  and  X_^  has  a Poisson  distribution  with 

parameter  X^  , i = l,2,...,k  . The  verification  that  Assumption  A is 

satisfied  for  this  case  is  deferred  to  Appendix  1. 

To  find  0*(x)  for  the  case  k=l  (where  all  the  subsystem  sample 

sizes  are  equal)  we  use  the  natural  ordering  of  x = x = 0,1,2,. . . and 

observe  that  0(A)  = X so  that  0*(x)  is  simply  the  value  of  X±  for 

which  P.  {X  < x}  = a . This  corresponds  to  the  standard  UMP  one-sided 
A1  ~ 

test  concerning  X^  , and  0*  coincides  with  the  "classical"  1-a  level 
confidence  bound  for  X^.  Tables  for  this  bound  have  usually  been  based 
on  the  relationship  between  the  Poisson  and  the  chi-squared  distributions. 

We  have  used  the  Poisson  distribution  directly  to  compute  Tables  1 and  2, 
which  give  0*(x)  for  a = .10  and  a = .05  for  values  of  x between 
0 and  149  . 
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For  the  case  k _>  2 ( k distinct  subsystem  sample  sizes)  the  deter- 
mination of  0*  requires  the  simultaneous  solution  of  k non-linear 
equations.  To  see  this  we  note  that  for  general  k 


(6) 


where 


= l px{x  = y>  > 

— - - 


(T) 


k 1 k y. 

P^{X=y}  = exp  X±|  ^ (A/Vy.l) 


i=l 


i=l 


JIV 

We  must  maximize  (6)  under  the  constraint  £ a.A.  = t and  obtain  0*(x) 

i=l  1 1 

as  the  value  of  t for  which  the  maximum  assumes  the  value  a . The 
constrained  maximization  may  be  accomplished  by  setting 


in  (6)  and  solving  simultaneously-  the  k-1  equations  obtained  by  setting 

the  partial  derivatives,  of  (6)  with,  respect  to  X^  , i = 1,2,..., k-1  , 

equal  to  zero,  with  the  kth  equation  obtained  by  setting  (6)  equal  to  a . 

This  is  equivalent  to  finding  the  solution  X*  = (X*  ,X*  , . . . ,A* ) with 

-L 

A*  > 0 , i = l,2,...,k  , of  the  simultaneous  equations 

1X7  pA{x ^ f + 1X7  px*x  — — * aXT  (*  -T  hh)  * 0 • J = 

J - k — 5 k \ i=l 

or 


(8) 


PA(X  < x)  - = 0 , J = 1,2,..., k-1  , 

0 — k dAk  — 
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and 


(9)  P-Jx  <x}  = a , 


and  setting 

(10)  0*tx)  = f a. A*  . 

“ A 1 1 


Observing  that 

(11)  3TTPA{X  = I}  = ( xt  “1)  PA{^=I}  ’ 

J — ' j ' — 

we  may  write  (8)  and  (9)  more  explicitly  as 


Ci2) 

and 


(13) 


l PA{x  = y}  = a 
y<x  — 


j = 1,2,... ,k-l  , 


where  P^{X_=  y}  is  given  by  (.7 ) • 

In  Appendix  2 it  is  shown  that  for  large  y = max(A  ,A„,. . . ,A,  ) the 

-L  xL 

asymptotically  optimal  upper  confidence  bound  for  0 ( A_)  is  given  by 


(HO 


I 


i=l 


a.X.  + c 
11  a 
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where  ea  is  the  upper  (lOOa)tl:1  percentage  point  of  the  standard  normal 
distribution.  This  bound  is  based  on  the  asymptotic  normal  distribution 
of  the  maximum  likelihood  estimate  of  0(A)  . We  therefore  choose  to 
generate  the  ordering  of  the  sample  points  x by  means  of  the  function 
(J)(x)  = 0 (x)  . This  guarantees  that  (i)  the  bound  0*(x)  will  be  asymp- 
totically optimal  and  (ii)  in  the  absence  of  strong  prior  information 
concerning  relative  subsystem  reliabilities,  8*(x)  will  be  "reasonable" 
for  small  values  of  y (and  hence  for  small  X^s)  in  the  sense  that  the 
subsystems  are  treated  evenhandedly  according  to  the  principle  of  maximum 
likelihood.  It  should  be  born  in  mind  that  only  the  ordering  of  the  sample 
points  is  determined  by  <Kx)  and  that  any  monotone  function  of  <j)  would 
produce  the  same  ordering.  Small  sample  deficiencies  in  the  maximum 
likelihood  bound  will  not  therefore  be  reflected  in  the  performance  of 
6* (x ) . 

For  the  case  k = 2 (two  distinct  subsystem  sample  sizes;).  Tables  3 
and  H give  0*(x)  for  a = .10  and  a = .05  for  the  first  100  ordered 
values  of  the  sample  point  X = (x  ,X2)  for  six  representative  values  of 
al  (-ai  = •-L(“),  -25  s -30,  1/3,  . 4o).  Note  that  a^  is  necessarily  less 

than  1/2  since  a^  <_  a2  = 1 - a^.  and  the  case  a = = 1/2  reduces 

to  the  case  k=l  . 

Example  1.  Consider  a weapons  system  consisting  of  five  subsystems 
and  suppose  that  the  entire  system  is  subjected  to  10,000  trials  (duty 
cycles)  with  necessary  repairs  being  made  after  each  subsystem  failure. 
Suppose  further  that  subsystem  number  five  is  redesigned  after  5,000  trials 
so  that  only  data  for  the  last  5,000  trials  of  that  subsystem  are  relevant 
to  the  reliability  of  the  final  system  configuration.  Suppose  that  the 
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observed  data  are  as  follows: 


Subsystem 

1 

2 

3 

It 

5 


Number  of 
Trials 

10,000 

10,000 

10,000 

10,000 

5,000 


Observed 

Failures 

1 

10 

5 

0 

1 


Since  the  first  four  subsystems  were  subjected  to  the  same  number  of 
trials  (n^  = 10,000)  we  combine  their  failure  data,  obtaining  X^  = l6  . 
For  the  remaining  subsystem  we  have  n^  = 5,000  and  X^  = 1 . To  find  a 
90%  lower  confidence  bound  for  MTBF  for  this  system  we  first  note  that 
a = 1/[10, 000(1/10, 000  + 1/5,000)]  = 1/3  by  (2).  From  Table  3 we  find 
that  the  upper  bound  .0*  for  0 corresponding  to  a^  = 1/3  and 
(X^jX^)  = (l6,l)  is  0*  = 8.422  . By  (3)  we  may  then  assert  that  at  the 
90%  level  of  confidence 

MTBF  _>  1/ [8. 422. Cl/10, 000.  + 1/5,000)]  = 396  . 

Example  2.  For  the  system  considered  in  Example  1,  suppose  that  the 
failure  data  were  the  same  as  those  given  above  except  that  the  redesign 
of  the  fifth  subsystem  was  not  performed  so  that  its  sample  size  is  also 
10,000  trials.  Since  all  the  sample  sizes  are  equal  we  are  in  the  case 
k=l  with  the  total  number  of  failures  being  X = IT  . For  a 90%  confidence 
bound  we  note  that  from  Table  1,  for  X = 17  the  value  of  the  upper  bound 
0*  is  23.606  . The  90%  confidence  statement  for  MTBF  then  becomes 

MTBF  >_  l/{ 23. 606(.1/10,000)J  = 424  . 
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It  should  be  noted  that  had  the  number  of  triala  in  these  examples 


been  aa  small  as  500  to  1000,  the  corresponding  analysis  would  still  be 
valid  since  the  Poisson  approximations  for  the  binomial  distributions 


would  still  be  very  close. 


3 . Approximations . In  Appendix  2 the  asymptotic  theory  (for  large 
Y = max(.A  ,X2, . ...  ,A  ) ) of  the  maximum  likelihood  ratio  procedure  for 
testing  hypotheses  concerning  0 is  used  to  obtain  approximate  1 - a 
level  confidence  bounds  for  0(A)  . Specifically,  if  jj  is  the  positive 
real  root  less  than  l/max ( ax , ag , . . . ,a^.)  of  the  equation 


(15) 


X 


1,2a 


= 2 


k 

l 

i=l 


ij + log(1"ai^ 


1 1— 


where  Xlj2a  is  the  upper  100 (2a)  percentage  point  of  the  chi-squared 
distribution  with  one  degree  of  freedom,  then  the  quantity 


(1 6) 


e(x) 


k 

l 


i=l 


a.X./ (l— a.p ) 
r r l 


is  an  approximate  upper  1 - a level  confidence  bound  for  0(A)  . This 
bound  is  asymptotically  optimal  and  asymptotically  equivalent  to  0 given 
by  (lU),  but  it  is  more  precise  than  the  latter  for  moderate  values  of  the 
s • Equation  (.15 ) may  easily  be  solved  for  p by  trial  and  error 
approximation  using  a hand  calculator.  For  the  case  k = 2 , 0 is  a 

good  approximation  to  6*  for  values  of  (X  ,X2)  beyond  those  tabulated. 
Comparisons  of  the  M.L.R.  approximations  with  the  exact  values  of  0*  from 
Tables  3 and  4 for  the  "largest"  pairs  (X_^,X2)  (corresponding  to  I = 100) 
for  each  case  are  given  in  Tables  A and  B below: 
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The  M.L.R.  hound  always  tends  to  underestimate  the  true  hound,  hut  the 
approximations  are  seen  to  he  within  10%  of  the  exact  values  in  all  cases 
and  within  5%  in  most  cases.  The  precision  of  0 will,  of  course,  improve 
further  for  larger  ' s . 

For  k >_  3 , if  some  of  the  X^'s  are  moderately  large,  it  is  reason- 
able to  expect  that  the  M.L.R.  hound  0 will  give  satisfactory  results. 

For  k >_  3 and  small  X^'s  a better  approximation  is  required.  We  observe 
that  the  asymptotic  hound  0 given  by  (l4)  depends  only  on  the  maximum 
likelihood  estimate  of  0 , i.e„, 

k 

(17)  9=1  a.X.  , 

u 11 
i=l 

and  the  maximum  likelihood  estimate  of  the  variance  of  0 , i.e., 

(18) 

i=. 


v 


- i 4* 


i i 


This  suggests  that  cases  involving  k 3 may  he  reduced  approximately  to 
corresponding  cases  with  k = 2 by  equating  the  corresponding  values  of 
0 and  v . A specific  procedure  for  accomplishing  this  may  he  described 
as  follows : 


(i)  Choose  an  integer  ^ , 1 <_  k < k , so  that  the  two  groups  of 

constants  a ,a  ,...,a  and  a 

-L  c.  & 

as  possible. 

(ii)  Compute  the  group  averages 


1+1’  Vj+2' 


a are  as  homogeneous 


a,  = 


_1 

k. 


*1 

l 

1 i=l 


a. 


k 

l *t  ■ 

i=kx+l  1 
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(iii)  Determine  the  equivalent  observations  X^  and  X^  for  the 
k = 2 case  by  setting 

aiXl  + a2X2  = 0 and  + S2X2  = ^ ’ 

where  0 and  v are  given  by  (17 ) and  (l8),  thereby  obtaining 

_ aT  6 - v _ v - a 0" 

Xi  = and  x2  = — — ± . 

al  ''a2-ai  ^ ai  ^a2~ai ) 

(iv)  Determine,  the  normalized  constants:  a*  , a*  for  the 

approximately  equivalent  k = 2 case  by 

a*  = a^/Ca^+aTg)  , a*  = a^/  (a^+a^ ) . 

(v)  Find  the  upper  1 - a confidence  bound  0*  for  k = 2 using 

al  ’ a2  ’ X1  and  X2  and  ^ or  Table  4,  interpolating  as  necessary. 

(vi)  Compute  the  1 - a upper  confidence  bound  0**  , say,  for  the 
original  (k  > 3)  problem  by  the  formula 

0**  = (J1+a2)e*  , 

where  0*  is  the  value  found  in  (v)  . 

The  validity  of  this  procedure  was  checked  for  a number  of  examples 
involving  k = 3 or  k = 4 using  the  M.L.R.  bound  9 computed  for  the 
original  example  and  the  corresponding  approximately  equivalent  k = 2 case 
(using  a^  and  instead  of  a*  and  a*  for  comparability).  The 
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agreement  in  the  corresponding  values  of  0 was  found  to  be  extraordinarily 
close,  usually  to  within  at  least  three  significant  figures.  The  procedure 
was  also  found  to  be  very  insensitive  to  the  choice  of  k . It  should  be 
noted  that  since  the  asymptotic  bound  depends  on  two  quantities  (F  and  v) 
it  is  not  reasonable  to  attempt  to  further  reduce  cases  involving  k > 2 
to  the  k = 1 case. 

The  proposed  procedure  is  illustrated  in  the  following  example: 

• Suppose  k = 3,  a^=l/6,  a^  = 1/3  , a3  = > the  observa- 
tions are  = 2 , = 3 , X^  = 5 > and  we  seek  an  upper  90%  confidence 


bound  for  0 . Then  from  (l7 ) and  (l8)  we  have  F = 3.8333  and  v = I.6389 
Choosing  k±  = 2 we  obtain  ^ = .25  and  a2  = ag  = .50  . 

From  the  formulas  of  (iii)  we  obtain  X1  = 4. 4444  and  x”2  = 5.4444  , and 
from  (iv),  a.^  = 1/3  and  a2  = 2/3  . Since  X^  and  X2  are  not  integers 
we  must  interpolate  in  Table  3 (a  = .10)  . For  ax  = 1/3  the  relevant 
entries  from  Table  3 are  as  follows: 


I 

fl 

Sl 

6* 

63 

4 

5 

7.339 

81 

4 

6 

8.176 

71 

5 

5 

7.592 

90 

5 

6 

8.533 

Four  point  linear  interpolation  yields 

9*  = (.5556)2(7.339)  + (.5556)(. 4444 ) (8. 176 + 7-592)  +( ,4444)2(8.533)  = 7.844 

Thus  from  (vi ) the  90%  upper  confidence  bound  for  0 is  given  by 

6**  = (7 .844) (.25 + .50)  = 5.883  . 
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The  values  of  the  M.L.R.  bound  0 for  the  original  k = 3 problem  and 
the  reduced  k = 2 problem  are  respectively  5-736  and  5-737  • The 
very  close  agreement  between  these  values  (even  though  they  both  underestimate 
0** ) strongly  supports  the  validity  of  the  suggested  procedure. 
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Appendix  1 


This  appendix  contains  the  proofs  of  the  four  general  propositions 
of  Section  2 and  the  verification  that  Assumption  A is  satisfied  for  the 
Poisson  confidence  hound  problem. 

Proof  of  Proposition  1:  Consider  e such  that  x Yr  y . Then  for 

all  \ e A , 

px(x  ^ 1 P^X^y)  » 

and  hence  for  all  t 

sup  P,  {Xy^x}  > sup  P,{Xc=y}  . 

AeS Ct ) — XeS(t)  - ” 

Thus,  recalling  the  definition  of  0*(x)  and  letting  S*  = S(0*(y))  = 

(_A:  0(A_)  = 0*(y)}  , we  have 

sup  P,  {x  ^ x)  > a , 

AeS*  - ~ 

so  that  0*(x)  _>  0*(yj  by  the  monotonicity  part  of  Assumption  A,  which 
completes  the  proof. 

Proof  of  Proposition  2:  For  fixed  arbitrary  A_  e A , let  t = 0(A)  . 

Then  if  0**  represents  any  inverse  of  0*  taking  values  in  the  equivalence 
classes  of  generated  by  the  ordering  ( ^ ) , we  have 

Px{0(A)  >_  0*(x)}  = P (x*  0**(t)} 

< sup  P,{X^0**(t)} 

AeS(t)  - 

< h(t ) = a , 
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by  the  definition  of  8*  . Hence 

px(0(x)  i e*(x)}  >_  px{8(x)  < e*(x)}  > 1 - a , 

and  the  proof  is  complete. 

Proof  of  Proposition  3:  Proceeding  contrapositively  we  assume  that  there 

exists  a point  x'  e such  that  0(x'  ) < 0*(x_' ) . Let 

S'  = S(0(x*  ))  = {X_:  0(_X)  = 0(x')}  • Then  by  Assumption  A and  the  definition 
& 

of  0 we  have 

sup  P,{X  4 x'}  > a , 

XeS'  - 

or  since  0 is  monotone  in  ( X? ) , 

sup  P,{0(X)  < 0(x’ )}  > a . 

XeS’  - 

Thus,  there  exists  a X_'  e S'  such  that 

PA,{0CX)  < 0(X» )}  > a , 

which  contradicts  the  assumption  that  0 is  a 1-a  level  confidence  bound 
and  completes  the  proof. 

Proof  of  Proposition  U:  For  any  1-a  level  confidence  bound  0 let 
0(x)  = sup  0(y)  . Then  0 is  monotone  in  the  ordering  (^)  and  if  it  is 

also  a 1-a  level  confidence  bound.  Proposition  3 guarantees  that 


- 2Q 


0(x)  > 0*(x)  for  all  x e 3^  which  is  the  desired  result.  Suppose  contra- 
positively  that  0 is  not  a 1-a  level  confidence  bound.  Then  there 
exists  a A°  e A such  that 


P (0(X)  < 0(A°)}  > a . 
A_°  — — — 

But  0(x)  < 0(x)  for  all  x , so  that 

P {0(X)  1 0(1°)}  > a , 

o 


which  contradicts  the  assumption  that  0 is  a 1-a  level  confidence 
bound  and  completes  the  proof. 

We  demonstrate  that  Assumption  A is  satisfied  for  the  Poisson  confidence 

bound  problem  under  consideration  by  means  of  the  following  two  lemmas: 

For  each  t let  A.  be  a value  of  A_  for  which  the  supremum  of  Assumption 

k 


A is  attained,  i.e.,  for  which  (12)  is  satisfied. 

A 


Let  0(x)  = J a.X.  , 
— • 11 
i=l  x x 


and  let  9(x)  be  as  defined  in  (ll+). 


Lemma  1.  For  t > 0 


9t  A 


P,  {X  ^ < 0 if  and  only  if 


E.  {0(X)|9(X)  < 0(x)}  < t . 


Proof : Suppressing  the  dependence  on  t in  the  notation  for  the  components 

k 

of  A = {A  , . . . ,X  ) and  recalling  that  V a. A.  = t we  see  that 
~ t 12  k 1 1 


9A . 


9t 


“ - l/a^  , i = l,2,...,k,  which  when  combined  with  (ll)  yields 


(A1) 


*>  - l l ^(i1-  3 h • 

y^x  j=i  j / -t 


21  - 


Now  from  (12),  summing  on  j yields 


k-1  a /y  \ k=l 

(A2)  l l 2 \ - lj  ^ {^=^}  = I I 


yxx  j=l  a , 


y<x  j=l 


(X  = y> 


From  (Al)  and  (A2)  we  obtain 


(A3)  (X^x)  = 


a i 

ij=i  aj  / y** 


k 


- 1)  pxt{A=y> 


By  symmetry  (A3)  must  also  hold  with  a.,  y.  , X.  replacing  a , y , X 

ill  K.  K K. 

respectively  for  i = 1,2,..., k-1  . Hence  multiplying  each  of  these  equa- 
tions by  the  corresponding  X^  and  summing  we  obtain 


(Alt)  ( I A 


i,  5tP,lx4x}  = 

1=1  / -t 


_1 

2 


l 

i=l  a j yxx_ 


l (9(y)  - t)  P,  {X  =y}  , 
-t 


which  yields  the  desired  result  upon  observing  that 


Ea  {9(x)|x^y}  = (px  {X  l 0(y)  PA  {X=y}  , 

— t ' — t ' y<x  — t 


and  recalling  that  x ^ y 0(x_)  <_  0(y)  . Note  that  since  t > 0 , 

k 

X,  ^ 0 and  hence  £ X.  > 0 . 
t i=l  1 

Now  consider  a random  vector  Z_ = (Z  jZ^, . . . ,Zn)  where  the  Z^'s  are 
independent  non-negative  random  variables.  Let  C be  the  hypercube  in 
k-dimensional  space  defined  for  c £ (0,°°)  by 


— jz_  — ( z^ , z^ , • • • , ) : 0 z ^ ^ c ’ ^ l,2,...,k|  . 
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For  any  measurable  set  R in  the  non— negative  orthant  of  k-space,  for 
v >_  0 let 

R.(v)  = R 0 {z:  z.=v} 

1 — 1 

for  i = 1,2, . . . ,k  , 

Lemma  2.1  If  R C C , P{Z_eR}  > 0 , and  if  v > v'  implies  R^(v)  R^Cv') 
for  i = l,2,...,k  , then 

E{Z . I Z e R>  < E{Z  Jz.  < c} 
i'—  — x.1  x— 

for  i = 1,2, ... ,k  . 

Proof:  Let  f.  (z)  , i = l,2,...,k  , be  the  probability  densities  of  the 

Z^ ' s with  respect  to  a dominating  measure  y = y^  x y2  x • » • x y^  . Then 

taking  i = 1 and  letting  y^k  = y2  x x y^  , u = Cu2,u^, . . . ,u^)  , 

P (c)  = P{Z  < c}  , and  f ^ ^(u)  = .II  f . (u.  ) , we  have 

(A5)  E{Z  |ZeR}  = pv/.yr  v f (v)  f (k_1 ) (u ) dy  ^ 5 dy  . 

1-  PiZeRi  jQ  1 JRi(v)  - 1 

But  the  inner  integral  on  the  right  is  descreasing  in  v so  that  by  the 
well  known  inequality  which  states  that  if  and  gg  are  two  functions 

monotone  in  opposite  senses  and  Y is  a random  variable,  then 
E g1(Y)  gg ( Y ) <_  E gx(Y)  E gg(Y)  , we  have 

1The  author  benefited  from  helpful  discussions  with  Rupert  G.  Miller 

concerning  the  proof  of  this  lemma. 
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(A6) 


R1(v) 


c 


, (k-1 ) / \ , (k-l)\  fl^  , 

^ iv  ) TJT)  ivi 


tc 


V 


P.Cc)  dyl 


o 'JVv) 


Jk-l),  ^ , (k-1 )\  fl^  , 

f (-}  dy  ) i^I)  dyi 


= E{Z1|zi£c}  P(ZeR}/p  (c)  . 


Hence  by  (A5)  and  (A 6) 


E{Z  |z_eR>  <_E{Z1|Z1<_c}  , 


and  the  proof  is  complete. 

We  observe  that  for  any  i , E{Z^|z^<_c}  < E Z.  , provided  that 

P{Z^<_c}  < 1 , which  implies  that  the  conditional  c.d.f.  of  Z^  , given 
{Z±  <_  c}  , is  strictly  greater  than  the  c.d.f.  of  Z^  on  the  interval 
( 0 , c ) . 

To  verify  that  Assumption  A is  satisfied  for  the  Poisson  model  under 

A 

consideration,  let  Z^  = ajt^  , i = l,2,...,k  , so  that  0 given  by  (l4) 
becomes 


0Cx) 


Let  z_  = . . . ,z^)  and  for  fixed  arbitrary  x.  e 36.  » let 

R = { z_ : z_^0  , g (z_)  £ 0 (x ) } . 


Since  ^ — - g (z_)  > 0 for  all  z__>  0 , we  conclude  that  for  v > v*  , 

££  R^v)  implies  z_  e Ri(v')  for  i = l,2,...,k  . Also,  for  sufficiently 
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so  that  the  conditions  of 


large  c > 0 , RCC  and  P,  {Z  £ R>  > 0 

A — 

Lemma  2 are  satisfied.  Also,  if  > 0 , then  Pfz^.  <_  c}  < 1 and  by 

the  remark  following  Lemma  2,  ElZ  jz  <c)  < E Z.  . Hence  noting  that 
_ k xii 

0(X)  = £ Z.  we  have  for  A ^ 0 , 

i=l  1 “ 

^ k 

^{0^)  (9(x)^0U)}  = l ExfZijZ£R} 

k 

< l E.  Z.  = E,  0(x)  . 
i=l  . 1 1 L ~ 


_ k 

Thus,  recalling  that  E,  0(X.)  = £ a.A.  — t , we  have  by  Lemma  1, 

-t  i=l  1 1 


3t  hx^  at  pa  . - — <0  * 

"t 


and  therefore  hx(t)  is  decreasing  in  t for  all  x £ H . The  existence 
of  a smallest  value  of  t^Cx)  follows  from  the  continuity  and  form  of 
the  left-hand  members  of  equations  (12)  and  (l3)  and  Assumption  A is  therefore 
satisfied. 
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Appendix  2 


This  appendix  contains  the  detailed  development  of  the  approximations 

employed  in  Section  3.  These  approximations  are  based  on  the  asymptotic 

properties  of  the  maximum  likelihood  ratio  (M.L.R.)  test  statistic. 

As  before,  let  X = (x^, . . . ,Xfc)  and  A = (A^, . . . ^ where 

for  each  i , X is  a Poisson  random  variable  with  parameter  A and 

i 

the  X.'s  are  independent.  To  obtain  the  M.L.R.  test  of  H : 0=6  vs. 

k 0 0 

Hl:  0^0O  ’ we  let  A0  = : 60  = iI1  ai^i^  and  observe  that  the  likelihood 

function  is 

k 

L(A_,X)  = exp{  £ (X.  log  A.  - A.  - log  X.!)}  . 

i=l  1 1 1 1 

Since  the  unconstrained  maximum  likelihood  estimate  (m.L.E.)  of  A.  is 
, i = l,2,...,k  , we  have 

k 

max  log  L(A,X)  = J (X.  log  X±  - X.  - log  X± ! ) . 

To  obtain  the  M.L.E. 's  under  the  constraint  that  A_ e AQ  Ci.e.,  HQ  is 
true)  we  use  the  Lagrange  multiplier  method  and  observe  that  for 
j — 1,2,. .,,k 

3 k 

— (log  L(X,A)  + y(  l a. A.  - 0 )}  = 0 

j i=l  1 1 U 


if  and  only  if 


X. 

+ a . y — 1 . 

A, 
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Thus  the  restricted  M.L.E.!s  are 


A,  = X / ( 1-a  y ) , 


A /\  , 


for  i = l,2,...,k  , with  y = y(0n)  determined  by 


(AT) 


l aX/(l-o)  = 0 . 

i=l  u 


Thus,  if  S = max  L/max  L , the  M.L.R.  statistic  is 

Aq  A 


(A8) 


T = -2  log  S 

(i 


/X.  \ X. 

. log  X.  - X.  - X.  log f:-"---  ) + —— 
11  1 i yi— a^p J l-ai 


k 

= 2 ^ 

k v ( a y \ 

- 2 h i^jr + • 


or 


(A9) 


T = 20 Qy  +2  l X.  log(l-a.G ) . 
i=l 


If  T were  based  on  n i.i.d.  observations  we  would  expect  T to 

be  asymptotically  chi-squared  with  one  degree  of  freedom  as  n ->-  oo  . 

For  the  present  case  (since  k is  fixed)  a slightly  different  analysis  is 

needed.  Let  Z.  = A?^  2(X„-A . ) , i = l,2,...,k  . Then  the  Z.’s  are 
1111  773  x 

bounded  in  probability  and  if  A_.  00  , then  is  asymptotically  a 

1/2 

standard  normal  variable.  Row  X.  = A.  Z.  + A.  and  from  (AT ) 

li  li 


k a.  A.1//2  Z.  + a. A.  k 

l I a. A. 

i=l  1 - a.y  0 i=l  1 1 


and  hence 


(A10) 


k a. (A. 1/2  Z.  + a. A.y) 

l i- = 0 . 

u /s 

i=l  1 - a.y 

i 
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Let  Y = maxCV-^Ag, . . . ,A  ) and  consider  Y large.  Then  using  the  symbol 
0 tg(-Y))  "to  represent  a random  quantity  which  when  divided  by  g(y)  is 
bounded  in  probability  as  Y ^ , we  observe  that  (A10)  implies 

U = Op(Y  1//2)  so  that  (l-a^y)  = 1 + 0 (y  ) for  each  i , and 

hence 


(AH) 


k 

1 a-  A. 
i-1  1 1 


1/2 


Z. 

1 


k (1  + o (y-1/2)> 

2,  ^ 


l atA . 

i=l  1 1 


^ ^ 12  2 ^/p 

Now  log  (l-a.y)  = -a.y  - ^ + 0p(-Y  ) , so  that  (A9)  becomes 


k k / 

y 1 aiXi  - Z X a y + -a  y + 
i=l  i=l  V 1 

f - v a , 1/2  _ 1 -2  r 2,  \ 

-w.iaixi  \-iv  J ,aixi} 

„ i-l  i=l  j 


%(y-3/2>)  ■ 

+ 0p(Y"1/2)  • 


Recalling  (All)  we  have 


T - 2 «yi+op(Y-1X2))  l a‘X.  I ay.i  ♦ 0 (y-±,2i 


(A12) 


i=l 


= y2  l a2A.  + 0 (y~1/2)  ■ 

i-1  1 1 P 


1 ~2  | 2, 

i=l 


-1/2. 


Hence,  letting  e represent  a generic  random  quantity  approaching  zero 
in  probability  as  y -*  «>  , we  have  by  (All)  and  (A12), 
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Thus,  for  large  y , T is  approximately  the  square  of  a standard  normal 

variate  and  the  approximate  two-sided  test  of  HQ:  0=0^  at  the  a 

2 2 

significance  level  is  to  reject  if  T > xn  » where  Xn  is  the 

l , ot  i , Ou 

th 

upper  a quantile  of  the  chi-squared  distribution  with  one  degree  of 
freedom.  The  corresponding  asymptotic  two-sided  confidence  statement  for 
0 is 


i=l 


a.X . - 
i i 


Ca/2 


< 0 < 


k 

I 

i=l 


a.X. 

l i 


1 th 

with  probability  1 - a , where  ca/2  ^he  uPPer  (77  a)  quantile 
of  the  standard  normal  distribution.  The  corresponding  one-sided  confidence 

A A 

statement  at  level  1 - a is  0 <_  0 where  0 is  given  by  (lU).  To  obtain 
an  asymptotic  upper  (one-sided)  level  1 - a confidence  bound  for  0 using 


the  M.L.R. 
interval 


statistic 
(O,  1/max (a 


directly,  one  first  finds  the  root 
^ja^, . . . ,a^) ) of  the  equation 


Pq  lying  in  the 
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{Alb ) 


T = 


2 

Xl,2a  5 


where  T is  given  by  (a8).  Substituting  y^  for  y in  (AT ) then  yields 
0q  = 0 , the  required  upper  1 - a confidence  bound  for  0 . The 
negative  -root  of  (Al4)  will  yield  a lower  confidence  bound  for  0 . 
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Appendix  3 


Tables  3 and  4 for  the  case  k.  = 2 were  produced  by  solving  (.12) 
and  (l3)  simultaneously  for  A*  and  A*  by  numerical  methods  and 
computing  6*  from  ClO).  The  numerical  methods,  used  involved  a finitized 
analogue 'of  Newton’s  method  for  two  dimensions.  The  calculations  were 
done  in  part  on  a self-contained  mini-computer  system  (Wang  2200)  and  in 
part  on  an  IBM  360/60  system.  In  a small  proportion  of  the  cases  considered 
the  procedure  failed  to  converge  In  a reasonable  number  of  iterations  and 
the  tabulated  values  were  obtained  by  interpolation  according  to  the  values 

✓s 

of  9 given  by  (l4)  subject  to  the  requirement  that  0*  be  monotone  in 
the  given  ordering. 

An  examination  of  Tables  3 and  4 shows  that  often  the  same  or  very 
nearly  the  same  value  of  0 occurs  for  several  successive  sample  points. 
This  reflects  the  fact  that  at  certain  stages  an  additional  pair 
(xl>x2)  = x may  contribute  very  little  to  the  quantity  P^{x  ^ x}  in  the 
vicinity  of  the  solution  pair  (A^,A*)  = A* . In  higher  dimensions  (k  > 2) 
this  phenomenon  might  be  exploited  to  condense  the  corresponding  tables 
considerably. 

Tables  1 and  2 were  computed  by  directly  generating  the  required 
Poisson  probabilities  rather  than  by  the  traditional  method  utilizing  the 
corresponding  chi-squared  probabilities. 

The  author  wishes  to  acknowledge  the  assistance  rendered  by  Mr.  David 
Pasta  in  the  performance  of  the  computations  described  above. 
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TABLE  1 


ALPHA- .10  K~1 


X 

THETA* 

X 

THETA* 

X 

THETA* 

0 

2.303 

50 

60.339 

100 

114.075 

1 

3.890 

51 

61.429 

101 

115.138 

O 

5.322 

52 

62.518 

102 

116.202 

3 

6.681 

53 

63.605 

103 

117.265 

4 

7.994 

54 

64.692 

104 

118.327 

5 

9.275 

55 

65.779 

105 

119.390 

6 

10.532 

56 

66.864 

106 

120.452 

7 

11.771 

57 

67.949 

107 

121.514 

8 

12.995 

58 

69.033 

108 

122.576 

9 

14.206 

59 

70.116 

109 

123.637 

10 

15.407 

60 

71,199 

110 

124.698 

11 

16.598 

61 

72.281 

111 

125.759 

12 

17.782 

62 

73.362 

112 

126.819 

13 

18.958 

63 

74.442 

113 

127.879 

14 

20.128 

64 

75.523 

114 

128.939 

15 

21.292 

65 

76.602 

115 

129.999 

16 

22.452 

66 

77.680 

116 

131.059 

17 

23.606 

67 

78.759 

117 

132.118 

18 

24.756 

68 

79.836 

118 

133.177 

19 

25,902 

69 

80.913 

119 

134.235 

20 

27.045 

70 

81.990 

120 

135.294 

21 

28.184 

71 

83.066 

121 

136.352 

22 

29.320 

72 

84.141 

122 

137.410 

23 

30.453 

73 

85.216 

123 

138.468 

24 

31,583 

74 

86.291 

124 

139.525 

25 

32.711 

75 

87.364 

125 

140.582 

26 

33 . 836 

76 

88.438 

126 

141.640 

27 

34,959 

77 

89.511 

127 

142,696 

28 

36.080 

78 

90.583 

128 

143.753 

29 

37.198 

79 

91.655 

129 

144.809 

30 

38,315 

80 

92.727 

130 

145.865 

31 

39.430 

81 

93.798 

131 

146.921 

32 

40,543 

82 

94.869 

132  . 

147.977 

33 

41.654 

83 

95.939 

133 

149.033 

34 

42.764 

84 

97.009 

134 

150.088 

35 

43,872 

85 

98.078 

135 

151,143 

36 

44.978 

86 

99,147 

136 

152.198 

37 

46.083 

87 

100.216 

137 

153.253 

38 

47.187 

88 

101.284 

138 

154.307 

39 

48.289 

89 

102.352 

139 

155.361 

40 

49.390 

90 

103.419 

140 

156.416 

41 

50,490 

91 

104.486 

141 

157.470 

42 

51.589 

92 

105.553 

142 

158.523 

43 

52.686 

93 

106.620 

143 

159.577 

44 

53.782 

94 

107.686 

144 

160.630 

45 

54.878 

95 

108.751 

145 

161,683 

46 

55.972 

96 

109.817 

146 

162.736 

47 

57.065 

97 

110.882 

147 

163.789 

48 

58.158 

98 

111.946 

148 

164.842 

49 

59.249 

99 

113,011 

149 

165.894 

i 


TABLE  2 


ALPHA" . 05 

X 

THETA* 

0 

2.996 

1 

4.744 

2 

6 . 296 

3 

7 . 754 

4 

9.153 

5 

10.513 

6 

11.842 

7 

13.148 

8 

14.434 

9 

15.705 

10 

. 16.962 

11 

18.207 

12 

19.443 

13 

20.669 

14 

21.886 

15 

23 . 097 

16 

24.301 

17 

25.499 

18 

26.692 

19 

27.879 

20 

29.062 

21 

30.240 

22 

31.415 

23 

32.585 

24 

33.752 

25 

34.916 

26 

36.076 

27 

37.234 

28 

38.389 

29 

39.541 

30 

40.690 

31 

41.837 

32 

42.982 

33 

44.125 

34 

45.266 

35 

46,404 

36 

47.541 

37 

48.675 

38 

49.808 

39 

50.940 

40 

52.069 

41 

53,197 

42 

54.324 

43 

55,449 

44 

56.572 

45 

57.695 

46 

58.816 

47 

59,935 

48 

61.054 

49 

62.171 

K“  1 


X 

THETA* 

50 

63.287 

51 

64.402 

52 

65.516 

53 

66.628 

54 

67.740 

55 

68,850 

56 

69.960 

57 

71.069 

58 

72.176 

59 

73.283 

60 

74.389 

61 

75.494 

62 

76.599 

63 

77.702 

64 

78.805 

65 

79.906 

66 

81.007 

67 

82.108 

68 

83.207 

69 

84.306 

70 

85,404 

71 

86,501 

72 

87.598 

73 

88.694 

74 

89,790 

75 

90.884 

76 

91.979 

77 

93.072 

78 

94,165 

79 

95.257 

80 

96.349 

81 

97,441 

82 

98.531 

83 

99.621 

84 

100,711 

85 

101.800 

86 

102,888 

87 

103.977 

88 

105.064 

89 

106.151 

90 

107,238 

91 

108.324 

92 

109.409 

93 

110.494 

94 

111.579 

95 

112.663 

96 

113.747 

97 

114.831 

98 

115.914 

99 

116,996 

X 

THETA* 

100 

118.078 

101 

119.160 

102 

120.241 

103 

121.322 

104 

122.403 

105 

123.484 

106 

124.564 

107 

125.643 

108 

126.722 

109 

127.801 

110 

128.879 

111 

129.957 

112 

131.035 

113 

132.112 

114 

133.189 

115 

134.265 

116 

135.342 

117 

136,418 

118 

137,493 

119 

138.568 

120 

139.643 

121 

140.718 

122 

141.792 

123 

142.866 

124 

143.940 

125 

145.014 

126 

146.087 

127 

147.160 

128 

148.232 

129 

149.305 

130 

150.376 

131 

151.449 

132 

152.520 

133 

153,592 

134 

154.663 

135 

155.733 

136 

156.804 

137 

157.873 

138 

158,944 

139 

160.014 

140 

161,083 

141 

162.152 

142 

163.221 

143 

164.290 

144 

165.359 

145 

166,427 

146 

167.495 

147 

168.563 

148 

169.630 

149 

170.698 

ii 


TABLE  3 


ALPHA- 

*10 

Al^.10 

I 

XI 

X2 

THETA* 

1 

0 

0 

2.072 

2 

1 

0 

2.079 

3 

2 

0 

2.111 

4 

3 

0 

2.163 

5 

4 

0 

2,225 

6 

5 

0 

2,293 

7 

6 

0 

2.367 

8 

7 

0 

2.444 

9 

8 

0 

2.524 

10 

9 

0 

2.606 

11 

10 

0 

2.690 

12 

11 

0 

2.776 

13 

12 

0 

2.863 

14 

13 

0 

2.951 

15 

14 

0 

3.040 

16 

15 

0 

3,130 

17 

0 

1 

3.501 

18 

16 

0 

3.501 

19 

1 

1 

3.507 

20 

17 

0 

3.507 

21 

Aw* 

l 

3,540 

w 

Am  A— 

18 

0 

3.540 

23 

3 

1 

3.592 

24 

19 

0 

3.592 

25 

4 

1 

3.654 

26 

20 

0 

3.654 

27 

5 

1 

3.723 

28 

21 

0 

3.723 

29 

6 

1 

3.797 

30 

22 

0 

3.797 

31 

7 

1 

3.875 

32 

8 

1 

3,955 

33 

23 

0 

3.955 

34 

9 

1 

4.038 

35 

24 

0 

4.038 

36 

10 

1 

4.123 

37 

25 

0 

4,123 

38 

11 

1 

4.209 

39 

26 

0 

4.209 

40 

12 

1 

4,297 

41 

27 

0 

4,297 

42 

0 

2 

4.790 

43 

13 

1 

4,790 

44 

28 

0 

4.790 

45 

1 

2 

4.796 

46 

14 

1 

4,796 

47 

29 

0 

4,796 

48 

2 

2 

4.829 

49 

15 

1 

4.829 

50 

30 

0 

4,829 

I 

XI 

X2 

THETA* 

51 

3 

2 

4.881 

52 

16 

1 

4,881 

53 

31 

0 

4.881 

54 

4 

Am 

4.943 

55 

17 

1 

4.943 

56 

32 

0 

4.943 

57 

5 

2 

5.013 

58 

18 

1 

5.013 

59 

33 

0 

5.013 

60 

6 

2 

5.087 

61 

19 

1 

5.087 

62 

34 

0 

5.087 

63 

7 

*5 

Am 

5.165 

64 

20 

1 

5,165 

65. 

35 

0 

5,165 

66 

8 

2 

5.246 

67 

21 

1 

5.246 

68 

36 

0 

5.246 

69 

9 

2 

5.329 

70 

22 

1 

5.329 

71 

37 

0 

5.329 

72 

10 

D 

Am 

5,414 

73 

23 

l 

5.414 

74 

11 

Am 

5.501 

75 

38 

0 

5.501 

76 

24 

1 

5.501 

77 

12 

o 

Am 

5,589 

78 

0 

3 

6.013 

79 

39 

0 

6.013 

80 

25 

1 

6.013 

81 

13 

Am 

6.013 

82 

1 

3 

6.019 

83 

40 

0 

6.019 

84 

26 

1 

6,019 

85 

14 

2 

6.019 

86 

■j 

3 

6.052 

87 

41 

0 

6,052 

88 

27 

1 

6.052 

89 

15 

o 

Am 

6.052 

90 

3 

3 

6,104 

91 

42 

0 

6.104 

92 

28 

1 

6.104 

93 

16 

2 

6.104 

94 

4 

3 

6,1 66 

95 

43 

0 

6.166 

96 

29 

1 

6 , 166 

97 

17 

2 

6,166 

98 

5 

3 

6.236 

99 

44 

0 

6.236 

100 

30 

1 

6,236 

iii 


TABLE  3 < CONT « ) 


ALPHAS 

.10 

Al= . 20 

I 

XI 

X2 

THETA* 

1 

0 

0 

1.842 

2 

1 

0 

1.872 

3 

o 

Am 

0 

1.979 

4 

3 

0 

2.123 

5 

4 

0 

2.286 

6 

5 

0 

2.460 

7 

0 

1 

3.111 

8 

6 

0 

3.111 

9 

1 

1 

3.142 

10 

7 

0 

3.142 

11 

n 

Am 

1 

3.251 

12 

8 

0 

3.251 

13 

3 

1 

3.398 

14 

9 

0 

3.399 

15 

4 

1 

3.567 

16 

10 

0 

3.567 

17 

5 

1 

3 . 748 

18 

11 

0 

3.751 

19 

0 

n 

Am 

4.258 

20 

6 

1 

4.258 

21 

1 

2 

4.288 

Am  4. 

12 

0 

4.288 

23 

7 

1 

4.288 

24 

2 

2 

4.398 

25 

13 

0 

4.398 

26 

8 

1 

4.398 

27 

3 

2 

4.548 

28 

14 

0 

4.548 

29 

9 

1 

4.548 

30 

4 

•> 

Am 

4.719 

31 

15 

0 

4.719 

32 

10 

1 

4.721 

33 

5 

2 

4.904 

34 

0 

3 

5.344 

35 

16 

0 

5.344 

36 

11 

1 

5.344 

37 

6 

2 

5.344 

38 

1 

3 

5.362 

39 

17 

0 

5.438 

40 

12 

1 

5.486 

41 

7 

o 

4. 

5.486 

42 

2 

3 

5.486 

43 

18 

0 

5.486 

44 

13 

1 

5.486 

45 

8 

2 

5.486 

46 

3 

3 

5.636 

47 

19 

0 

5.636 

48 

14 

1 

5.636 

49 

9 

2 

5.637 

50 

4 

3 

5.810 

I 

XI 

X2 

THETA* 

51 

20 

0 

5.810 

52 

15 

1 

5.810 

53 

0 

4 

6 . 395 

54 

10 

2 

6.395 

55 

5 

3 

6.395 

56 

21 

0 

6.395 

57 

16 

1 

6 . 395 

58 

. 1 

4 

6.414 

59 

11 

2 

6.431 

60 

6 

3 

6.434 

61 

22 

0 

6.536 

62 

17 

1 

6.536 

63 

9 

4 

6.536 

64 

12 

2 

6.536 

65 

7 

3 

6.536 

66 

23 

0 

6.536 

67 

18 

1 

6.536 

68 

3 

4 

6.686 

69 

8 

3 

6.688 

70 

13 

2 

6.688 

71 

24 

0 

6,688 

72 

4 

4 

6.857 

73 

19 

1 

6.857 

74 

9 

3 

6.864 

75 

14 

o 

Am 

6.867 

76 

25 

0 

7.027 

77 

0 

5 

7.420 

78 

5 

4 

7.420 

79 

20 

1 

7.420 

80 

10 

3 

7.420 

81 

15 

*3 

Am 

7.420 

82 

26 

0 

7,420 

83 

1 

5 

7.420 

84 

6 

4 

7.463 

85 

21 

1 

7,479 

86 

11 

3 

7.490 

87 

16 

aL. 

7.498 

88 

2 

5 

7.561 

89 

27 

0 

7.561 

90 

7 

4 

7.562 

91 

22 

1 

7.562 

92 

12 

3 

7.562 

93 

17 

7.562 

94 

3 

5 

7.712 

95 

28 

0 

7.712 

96 

8 

4 

7.714 

97 

23 

1 

7.714 

98 

13 

3 

7.714 

99 

18 

2 

7.714 

100 

4 

5 

7.884 

iv 


TABLE  3 <CONT.  ) 


ALPHA=  .10  Al=.25  K=2 


I 

XI 

X2 

THETA# 

1 

0 

0 

1.727 

2 

1 

0 

1.781 

3 

2 

0 

1.944 

4 

3 

0 

2.151 

5 

4 

0 

2.381 

6 

0 

1 

2.917 

7 

5 

0 

2.917 

8 

1 

1 

2.972 

9 

6 

0 

2.972 

10 

o 

Am 

1 

3.142 

11 

7 

0 

3.142 

12 

3 

1 

3.362 

13 

0 

o 

Am 

3.991 

14 

4 

1 

3.991 

15 

8 

0 

3.991 

16 

1 

2 

4.047 

17 

5 

1 

4.047 

18 

9 

0 

4.047 

19 

2 

2 

4.217 

20 

6 

1 

4.220 

21 

10 

0 

4.220 

oo 

Am  Am 

3 

2 

4.433 

23 

7 

1 

4.450 

24 

11 

0 

4.450 

25 

0 

3 

5.011 

26 

4 

2 

5.011 

27 

8 

1 

5.011 

28 

12 

0 

5.011 

29 

1 

3 

5.066 

30 

5 

2 

5.066 

31 

9 

1 

5.066 

32 

13 

0 

5.066 

33 

2 

3 

5.238 

34 

6 

o 

Am 

5.242 

35 

10 

1 

5.242 

36 

14 

0 

5.242 

37 

3 

3 

5.456 

38 

7 

2 

5.489 

39 

0 

4 

5.995 

40 

11 

1 

5.995 

41 

15 

0 

5.995 

42 

4 

3 

5.995 

43 

8 

n 

Am 

5.995 

44 

1 

4 

6.051 

45 

12 

1 

6.051 

46 

16 

0 

6.051 

47 

5 

3 

6.051 

48 

9 

2 

6.051 

49 

2 

4 

6.224 

50 

13 

1 

6.224 

I 

XI 

X2 

THETA# 

51 

17 

0 

6 . 224 

52 

6 

3 

6.230 

53 

10 

o 

Am 

6.230 

54 

3 

4 

6.443 

55 

14 

1 

6.443 

56 

18 

0 

6.443 

57 

7 

3 

6.498 

58 

0 

5 

6.956 

59 

11 

2 

6.956 

60 

4 

4 

6.956 

61 

15 

1 

6.956 

62 

8 

3 

6.956 

63 

19 

0 

6.956 

64 

1 

5 

6.986 

65 

12 

2 

7.028 

66 

5 

4 

7.028 

67 

16 

1 

7.028 

68 

9 

3 

7.028 

69 

20 

0 

7.028 

70 

o 

Am 

5 

7.044 

71 

13 

o 

Am 

7.139 

72 

6 

4 

7.215 

73 

17 

1 

7.215 

74 

10 

3 

7.286 

75 

21 

0 

7.294 

76 

3 

5 

7.407 

77 

14 

o 

Aw 

7.407 

78 

7 

4 

7.449 

79 

0 

6 

7.899 

80 

18 

1 

7.899 

81 

11 

3 

7.899 

82 

4 

5 

7.899 

83 

on 

A M Am 

0 

7.899 

84 

15 

2 

7.899 

85 

8 

4 

7.899 

86 

1 

6 

7.900 

87 

19 

1 

7.914 

88 

12 

3 

7.914 

89 

5 

5 

7.914 

90 

23 

0 

7.914 

91 

16 

o 

Am 

7,914 

92 

2 

6 

8,130 

93 

9 

4 

8,130 

94 

20 

1 

8.130 

95 

13 

3 

8.130 

96 

6 

5 

8.136 

97 

24 

0 

8.136 

98 

17 

2 

8,136 

99 

3 

6 

8,351 

100 

10 

4 

8.355 

v 


TABLE  3 ( CONT. ) 


ALPHA-  .10  A1 = . 30  K=2 


I 

XI 

X2 

THETA* 

1 

0 

0 

1.612 

2 

1 

0 

1.703 

3 

2 

0 

1*941 

4 

3 

0 

2.229 

5 

0 

1 

2.723 

6 

4 

0 

2.723 

7 

1 

1 

2.819 

8 

rt 

Am 

1 

3.054 

9 

5 

0 

3.074 

10 

0 

2 

3.726 

11 

3 

1 

3.726 

12 

6 

0 

3.726 

13 

1 

2 

3.820 

14 

4 

1 

3.823 

15 

7 

0 

3.823 

16 

2 

■9 

4.059 

17 

5 

1 

4.094 

18 

8 

0 

4.096 

19 

0 

3 

4.676 

20 

3 

2 

4.676 

21 

6 

1 

4.676 

oo 

Am 

9 

0 

4,676 

23 

1 

3 

4.772 

24 

4 

9 

a-. 

4.776 

OCT 

w 

7 

1 

4.776 

26 

10 

0 

4.776 

27 

2 

3 

5,012 

28 

5 

2 

5.060 

29 

8 

1 

5.087 

30 

11 

0 

5.101 

31 

0 

4 

5.597 

32 

3 

3 

5.597 

33 

6 

9 

Am 

5.597 

34 

9 

1 

5.597 

35 

12 

0 

5.597 

36 

1 

4 

5.691 

37 

4 

3 

5,696 

38 

7 

9 

5.696 

39 

10 

1 

5,696 

40 

2 

4 

5.932 

41 

13 

0 

5.932 

42 

5 

3 

5.991 

43 

8 

9 

Am 

6.040 

44 

0 

5 

6.492 

45 

11 

1 

6.492 

46 

3 

4 

6.492 

47 

14 

0 

6.492 

48 

6 

3 

6.492 

49 

9 

2 

6.492 

50 

1 

5 

6.589 

I 

XI 

X2 

THETA* 

51 

12 

1 

6.589 

52 

4 

4 

6.594 

53 

15 

0 

6.594 

54 

7 

3 

6.594 

55 

10 

o 

Am 

6.594 

56 

A— 

5 

6 . 830 

57 

13 

1 

6.830 

58 

5 

4 

6.898 

59 

16 

0 

6.898 

60 

8 

3 

6.962 

61 

0 

6 

7.372 

62 

11 

2 

7.372 

63 

3 

5 

7,372 

64 

14 

1 

7.372 

65 

6 

4 

7.372 

66 

17 

0 

7.372 

67 

1 

6 

7.469 

68 

9 

3 

7.469 

69 

12 

2 

7.469 

70 

4 

5 

7.475 

71 

15 

1 

7.475 

72 

7 

4 

7.475 

73 

18 

0 

7.475 

74 

r> 

A*. 

6 

7.711 

75 

10 

3 

7.711 

76 

13 

2 

7.711 

77 

5 

5 

7.787 

78 

0 

7 

8.239 

79 

16 

1 

8.239 

80 

8 

4 

8.239 

81 

19 

0 

8.239 

82 

3 

6 

8.239 

83 

11 

3 

8.239 

84 

6 

5 

8.239 

85 

14 

2 

8.239 

86 

1 

7 

8.337 

87 

17 

1 

8,337 

88 

9 

4 

8,337 

89 

20 

0 

8.337 

90 

4 

6 

8.344 

91 

12 

3 

8,344 

92 

7 

5 

8.455 

93 

15 

2 

8,462 

94 

2 

7 

8.579 

95 

18 

1 

8.579 

96 

10 

4 

8,579 

97 

5 

6 

8 , 661 

98 

21 

0 

8.661 

99 

13 

3 

8.661 

100 

0 

8 

9.096 

vi 


TABLE  3 < CONT  * ) 


ALPHA= 

.10 

A1 =1/3 

I 

XI 

X2 

THETA# 

1 

0 

0 

1,534 

2 

1 

0 

1 .663 

3 

o 

0 

1.965 

4 

0 

1 

2.593 

5 

3 

0 

2.593 

6 

1 

1 

2.727 

7 

4 

0 

2.734 

8 

o 

1 

3.031 

9 

0 

2 

3.548 

10 

5 

0 

3.548 

11 

3 

1 

3.548 

12 

1 

2 

3.685 

13 

6 

0 

3.685 

14 

4 

1 

3.775 

15 

2 

2 

3.991 

16 

7 

0 

3.994 

17 

0 

3 

4.454 

18 

5 

1 

4 , 454 

19 

3 

2 

4.454 

20 

1 

3 

4.593 

21 

8 

0 

4.593 

*30 

**•  Am 

6 

1 

4.593 

23 

4 

2 

4.718 

24 

2 

3 

4.900 

25 

9 

0 

4,901 

26 

0 

4 

5.329 

27 

7 

1 

5.329 

28 

5 

2 

5.329 

29 

3 

3 

5.329 

30 

10 

0 

5.329 

31 

1 

4 

5.469 

32 

8 

1 

5.469 

33 

6 

*:> 

Am 

5.470 

34 

4 

3 

5.592 

35 

11 

0 

5,618 

36 

o 

Am 

4 

5 , 778 

37 

9 

1 

5,779 

38 

0 

5 

5,820 

39 

7 

o 

Am 

5.834 

40 

5 

3 

5.902 

41 

12 

0 

5,902 

42 

3 

4 

6.130 

43 

10 

1 

6.131 

44 

1 

5 

6.324 

45 

8 

2 

6.324 

46 

6 

3 

6.325 

47 

13 

0 

6.351 

48 

4 

4 

6,488 

49 

11 

1 

6.490 

50 

2 

5 

6.634 

I 

XI 

X2 

THETA* 

51 

0 

6 

6,668 

52 

9 

2 

6.668 

53 

7 

3 

6,668 

54 

14 

0 

6,720 

55 

5 

4 

6.755 

56 

12 

1 

6 . 755 

57 

3 

5 

6.984 

58 

1 

6 

7.164 

59 

10 

2 

7.164 

60 

8 

3 

7.164 

61 

15 

0 

7.164 

62 

6 

4 

7.164 

63 

4 

5 

7.339 

64 

13 

1 

7.339 

65 

o 

Am 

6 

7.474 

66 

11 

2 

7.474 

67 

0 

7 

7.515 

68 

9 

3 

7,515 

69 

16 

0 

7.515 

70 

7 

4 

7.515 

71 

5 

5 

7,592 

72 

14 

1 

7.592 

73 

3 

6 

7.671 

74 

12 

2 

7,702 

75 

1 

7 

7.990 

76 

10 

3 

7.990 

77 

8 

4 

7.990 

78 

17 

0 

7.990 

79 

6 

5 

7.991 

80 

15 

1 

7.991 

81 

4 

6 

8.176 

82 

13 

Am 

8.176 

83 

*3 

Am 

7 

8,302 

84 

11 

3 

8.302 

85 

0 

8 

8,310 

86 

9 

4 

8,310 

87 

18 

0 

8 . 320 

88 

7 

5 

8.347 

89 

16 

1 

8,422 

90 

5 

6 

8.533 

91 

14 

o 

8.533 

92 

3 

7 

8.648 

93 

12 

3 

8,648 

94 

1 

8 

8,807 

95 

10 

4 

8.807 

96 

19 

0 

8.807 

97 

8 

5 

8.807 

98 

17 

1 

8,807 

99 

6 

6 

8.808 

100 

15 

2 

8.808 

vii 


TABLE  3 (CONT.) 


ALPHA=  ,10  Al= . 40  K=2 


I 

XI 

X2 

THETA* 

1 

0 

0 

1.381 

2 

1 

0 

1.640 

3 

0 

1 

2.334 

4 

2 

0 

2.334 

5 

1 

1 

2.563 

6 

3 

0 

2.681 

7 

0 

o 

A. 

2.938 

8 

n 

Am 

1 

3.079 

9 

4 

0 

3.079 

10 

1 

2 

3.416 

11 

3 

1 

3.532 

12 

0 

3 

3.615 

13 

5 

0 

3.615 

14 

2 

2 

3.799 

15 

4 

1 

3.911 

16 

1 

3 

4 . 228 

17 

6 

0 

4.228 

18 

3 

2 

4.352 

19 

0 

4 

4.480 

20 

5 

1 

4.480 

21 

2 

3 

4.598 

22 

7 

0 

4.598 

23 

4 

A- 

4.814 

24 

1 

4 

5.014 

25 

6 

1 

5.014 

26 

3 

3 

5.144 

27 

8 

0 

5.144 

28 

0 

5 

5.221 

29 

5 

2 

5.239 

30 

2 

4 

5.477 

31 

7 

1 

5.477 

32 

4 

3 

5.598 

33 

9 

0 

5.599 

34 

1 

5 

5.781 

35 

6 

2 

5.782 

36 

'3 

4 

5.915 

37 

8 

1 

5.917 

38 

0 

6 

5 « 963 

39 

5 

3 

6.056 

40 

10 

0 

6.058 

41 

*3 

5 

6.238 

42 

7 

2 

6.238 

43 

4 

4 

6.364 

44 

9 

1 

6.366 

45 

1 

6 

6.535 

46 

6 

3 

6.535 

47 

11 

0 

6.535 

48 

3 

5 

6.672 

49 

8 

2 

6.674 

50 

0 

7 

6.698 

I 

XI 

X2 

THETA* 

51 

5 

4 

6.767 

52 

10 

1 

6.814 

53 

2 

6 

6.986 

54 

7 

3 

6.986 

55 

12 

0 

6.986 

56 

4 

5 

6.987 

57 

9 

n 

Am 

7.134 

58 

1 

7 

7.134 

59 

6 

4 

7.278 

60 

11 

1 

7.278 

61 

3 

6 

7.418 

62 

8 

3 

7.421 

63 

0 

3 

7.432 

64 

13 

0 

7.436 

65 

5 

5 

7.511 

66 

10 

2 

7.511 

67 

2 

7 

7.724 

68 

7 

4 

7.725 

69 

12 

1 

7.725 

70 

4 

6 

7.859 

71 

9 

3 

7.867 

72 

1 

8 

7.889 

73 

14 

0 

7.895 

74 

6 

5 

8.012 

75 

11 

o 

Am 

8.012 

76 

3 

7 

8.154 

77 

0 

9 

8.179 

78 

8 

4 

8.179 

79 

13 

1 

8.246 

80 

5 

6 

8.246 

81 

2 

8 

8.454 

82 

10 

3 

8.454 

83 

15 

0 

8.454 

84 

7 

5 

8.455 

85 

4 

7 

8.592 

86 

12 

2 

8.596 

87 

1 

9 

8.737 

88 

9 

4 

8.750 

89 

6 

6 

8.751 

90 

14 

1 

8.752 

91 

3 

8 

8.882 

92 

11 

3 

8.882 

93 

0 

10 

8.891 

94 

8 

5 

8.891 

95 

16 

0 

8.891 

96 

5 

7 

8.973 

97 

13 

2 

8.973 

98 

2 

9 

9.178 

99 

10 

4 

9.178 

100 

7 

6 

9.178 

viii 


TABLE  4 


ALPHA= 

* 05 

A1-.10 

r 

XI 

X2 

THETA* 

i 

0 

0 

2,695 

n 

1 

0 

2.700 

3 

A- 

0 

2.735 

4 

3 

0 

2.787 

5 

4 

0 

2.849 

6 

5 

0 

2.917 

7 

6 

0 

2.991 

8 

7 

0 

3,068 

9 

8 

0 

3.148 

10 

9 

0 

3.230 

11 

10 

0 

3.314 

12 

11 

0 

3.400 

13 

12 

0 

3.487 

14 

13 

0 

3.575 

15 

14 

0 

3.664 

16 

15 

0 

3.754 

17 

16 

0 

3.844 

18 

17 

0 

3.936 

19 

0 

1 

4.269 

20 

1 

1 

4.275 

21 

18 

0 

4.275 

oo 

Am  Am 

*Cm 

1 

4.309 

23 

19 

0 

4.309 

24 

3 

1 

4.360 

25 

20 

0 

4.360 

26 

4 

1 

4.422 

27 

21 

0 

4.422 

28 

5 

1 

4.491 

29 

'3'? 

Am  An 

0 

4.491 

30 

6 

1 

4,565 

31 

23 

0 

4.565 

32 

7 

1 

4.643 

33 

24 

0 

4.643 

34 

8 

1 

4,723 

35 

25 

0 

4.723 

36 

9 

1 

4.806 

37 

26 

0 

4,806 

38 

10 

1 

4.891 

39 

27 

0 

4.891 

40 

11 

1 

4.977 

41 

28 

0 

4.977 

42 

12 

1 

5.065 

43 

29 

0 

5.065 

44 

13 

1 

5,154 

45 

0 

2 

5.666 

46 

30 

0 

5,666 

47 

14 

1 

5.666 

48 

.1 

2 

5.672 

49 

15 

1 

5,672 

50 

31 

0 

5.672 

I 

XI 

X2 

THETA* 

51 

o 

Am 

2 

5.705 

52 

16 

1 

5.705 

53 

32 

0 

5.705 

54 

3 

2 

5.757 

55 

17 

1 

5.757 

56 

33 

0 

5.757 

57 

4 

2 

5,819 

58 

18 

1 

5.819 

59 

34 

0 

5.819 

60 

5 

2 

5.889 

61 

19 

1 

5.889 

62 

35 

0 

5.889 

63 

6 

o 

Am 

5.963 

64 

20 

1 

5.963 

65 

36 

0 

5.963 

66 

7 

2 

6.041 

67 

21 

1 

6.041 

68 

37 

0 

6.041 

69 

8 

o 

Am 

6.121 

70 

O'? 

Am  Am 

1 

6.121 

71 

38 

0 

6.121 

72 

9 

2 

6.204 

73 

23 

1 

6.204 

74 

39 

0 

6.204 

75 

10 

o 

Am 

6.289 

76 

24 

1 

6.289 

77 

40 

0 

6.289 

78 

11 

2 

6.376 

79 

25 

1 

6.376 

80 

41 

0 

6.376 

81 

12 

o 

Am 

6 * 464 

82 

26 

1 

6.464 

83 

0 

3 

6.978 

84 

42 

0 

6.978 

85 

13 

o 

Am 

6.978 

86 

27 

1 

6.978 

87 

1 

3 

6.984 

88 

43 

0 

6.984 

89 

14 

Am 

6.984 

90 

28 

1 

6.984 

91 

2 

3 

7,018 

92 

15 

*? 

Am 

7.018 

93 

44 

0 

7,018 

94 

29 

1 

7.018 

95 

3 

3 

7.069 

96 

16 

2 

7.069 

97 

45 

0 

7.069 

98 

30 

1 

7.069 

99 

4 

3 

7,132 

100 

17 

2 

7.132 

ix 


TABLE  4 ( CQNT « ) 


ALPHA= 

.05 

Ai= . 20 

I 

XI 

X2 

THETA* 

1 

0 

0 

2.396 

o 

1 

0 

2.425 

3 

2 

0 

2.534 

4 

3 

0 

2.677 

5 

4 

0 

2.840 

6 

5 

0 

3.015 

7 

6 

0 

3.198 

e 

0 

1 

3.795 

9 

7 

0 

3.795 

10 

1 

1 

3.825 

11 

8 

0 

3.825 

12 

r> 

Am 

1 

3.934 

13 

9 

0 

3.934 

14 

3 

1 

4.081 

15 

10 

0 

4.081 

16 

4 

1 

4.248 

17 

11 

0 

4.248 

18 

5 

1 

4.429 

19 

0 

2 

5.037 

20 

12 

0 

5.037 

21 

6 

1 

5.037 

22 

1 

2 

5.067 

23 

7 

1 

5.067 

24 

13 

0 

5.067 

25 

2 

2 

5.176 

26 

8 

1 

5.176 

27 

14 

0 

5.176 

28 

3 

2 

5.325 

29 

9 

1 

5.325 

30 

15 

0 

5.325 

31 

4 

2 

5.494 

32 

10 

1 

5.495 

33 

16 

0 

5 . 495 

34 

5 

2 

5.677 

35 

0 

3 

6.202 

36 

11 

1 

6.202 

37 

17 

0 

6.202 

38 

6 

2 

6.202 

39 

1 

3 

6.298 

40 

12 

1 

6.343 

41 

18 

0 

6.343 

42 

7 

2 

6.343 

43 

o 

Am 

3 

6.343 

44 

13 

1 

6.343 

45 

19 

0 

6.343 

46 

8 

*> 

6.343 

47 

3 

3 

6.493 

48 

14 

1 

6.493 

49 

20 

0 

6.493 

50 

9 

2 

6.493 

I 

XI 

X2 

THETA* 

51 

4 

3 

6 . 664 

52 

15 

1 

6 . 664 

53 

21 

0 

6 * 664 

54 

10 

o 

6 ♦ 666 

55 

5 

3 

6.850 

56 

0 

4 

7.323 

57 

16 

1 

7.323 

58 

22 

0 

7.323 

59 

11 

2 

7.323 

60 

6 

3 

7.323 

61 

1 

4 

7.365 

62 

17 

1 

7.411 

63 

23 

0 

7.464 

64 

12 

r> 

Am 

7.464 

65 

7 

3 

7,464 

66 

2 

4 

7.464 

67 

18 

1 

7.464 

68 

13 

r> 

Am 

7.464 

69 

24 

0 

7.464 

70 

8 

3 

7.464 

71 

3 

4 

7.614 

72 

19 

1 

7.614 

73 

14 

o 

Am 

7.614 

74 

25 

0 

7.614 

75 

9 

3 

7.614 

76 

4 

4 

7.787 

77 

20 

1 

7.787 

78 

15 

2 

7.787 

79 

26 

0 

7,787 

80 

10 

3 

7.790 

81 

5 

4 

7.975 

82 

0 

5 

8,410 

83 

21 

1 

8.410 

84 

16 

2 

8.410 

85 

27 

0 

8,410 

86 

11 

3 

8,410 

87 

6 

4 

8.410 

88 

1 

5 

8.420 

89 

22 

1 

8.499 

90 

17 

2 

8,552 

91 

28 

0 
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92 

12 

3 
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93 

7 

4 
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94 

2 
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23 
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8.552 

96 

18 
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97 

13 

3 

8.552 

98 

29 
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99 

8 

4 
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TABLE  4 < CONI . ) 


ALPHA-  « 05  Al=.25  K=2 
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XI 
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THETA* 

1 

0 

0 

2.247 

o 
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1 

0 

2.301 

3 

2 
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2.464 

4 

3 
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2.671 

5 
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0 
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3.557 

7 

5 
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3.557 
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9 
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10 

o 
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3.781 

11 
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3.781 
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3.997 

13 
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14 
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0 
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4.720 

17 

5 
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4.720 

18 

1 
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4.777 

19 

10 
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4.777 

20 

2 

o 
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4.946 

21 

6 

1 

4.948 
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A—  «L 

11 
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4.948 

23 

3 
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5.160 

24 

7 

1 

5.172 

25 

0 

3 

5.815 

26 

12 

0 

5.815 

27 

4 

2 

5.815 

28 

8 

1 

5.815 

29 

1 

3 

5.871 

30 

5 

2 

5.871 

31 

13 

0 

5.871 

32 

9 

1 

5.871 

33 
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3 

6.041 

34 

6 

o 
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6.044 

35 

14 

0 

6,044 

36 

10 

1 

6.044 

37 

3 

3 

6.257 

38 

7 

2 

6.274 

39 

11 

1 

6,274 

40 

15 

0 

6.274 

41 

0 

4 

6.864 

42 

4 

3 

6,864 

43 
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An 

6.864 

44 

12 

1 

6.864 

45 

16 

0 

6 . 864 

46 

1 
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6,921 

47 

5 
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6,921 

48 

9 
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49 

13 
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50 

17 

0 
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I 
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51 

2 
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52 

6 

3 

7.096 

53 

10 

2 
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54 

14 

1 
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55 

18 

0 
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56 

3 

4 

7.309 

57 

7 

3 

7.332 

58 

11 

*? 

Am 

7,409 

59 

0 

5 

7.884 

60 

15 

1 

7.884 

61 

19 

0 

7.884 

62 

4 

4 

7.884 

63 

8 

3 

7.884 

64 

12 

o 

Am 

7.884 

65 

1 

5 

7.941 

66 

16 

1 

7.941 

67 

20 

0 

7.941 

68 

5 

4 

7.941 

69 

9 

3 

7,941 

70 

13 

2 

7.941 

71 

2 

5 

8.113 

72 

17 

1 

8.113 

73 

21 

0 

8.113 

74 

6 

4 

8,118 

75 

10 

3 

8.140 

76 

14 

2 

8,175 

77 

3 

5 

8.175 

78 

18 

1 

8.175 

79 

22 

0 

8.221 

80 

7 

4 

8.367 

81 

0 

6 

8.882 

82 

11 

3 

8.882 

83 

15 

2 

8.882 

84 

4 

5 

8.882 

85 

19 

1 
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86 
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87 

8 
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1 
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8,938 

89 
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3 

8,938 

90 
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91 

5 

5 
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97 
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9,116 
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I 

XI 

X2 

THETA* 
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XI 
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THETA* 

1 

0 

0 
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0 
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53 

4 

4 

7.459 
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6 

4* 
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60 

8 

3 

7.777 
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1 
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8.289 

15 
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4.504 

65 
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0 

8.289 
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2 

o 

Am 

4.740 

66 

6 

4 

8.289 

17 

5 

1 

4.765 

67 

9 

3 

8.289 

18 

8 

0 

4.765 

68 

12 

2 

8.289 

19 

0 

3 

5.427 

69 

1 

6 

8.386 

20 

3 

2 

5.427 

70 

15 

1 

8.386 

21 

6 

1 

5.427 

71 

4 

5 

8.390 

90 

A*.  Am 
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5.427 

72 
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5.525 

74 
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25 
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5.525 
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26 
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76 
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27 

2 
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5.761 

77 

16 
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28 
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2 

5.796 

78 
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8.683 

29 

8 
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5.798 

79 

19 

0 
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30 

11 

0 

5.798 

80 

8 
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8,728 
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0 
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6.407 

81 

0 
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9,203 

32 

3 

3 

6.407 

82 

11 

3 

9,203 

33 
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Am 

6.407 

83 

3 

6 

9,203 

34 

9 

1 

6.407 

84 

14 

2 

9,203 

35 

12 

0 

6.407 

85 

17 

1 

9.203 

36 

1 

4 

6.503 

86 

6 

5 

9.203 

37 

4 

3 

6.506 

87 

20 

0 

9.203 

38 

7 

Am 

6.506 

88 

9 

4 

9.203 

39 

10 

1 

6.506 

89 

1 

7 

9.300 

40 

13 

0 

6.506 

90 

12 

3 
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2 

4 

6.742 

91 

4 
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9,305 

42 

5 

3 

6.785 

92 
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9.305 

43 
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Am 

6.800 

93 

7 

5 

9.305 

44 

11 

1 

6.809 

94 

18 

1 

9,305 

45 

14 

0 

6.810 

95 

21 

0 

9.305 

46 

0 

5 

7.358 

96 

10 

4 

9.305 

47 

3 

4 

7.358 

97 

2 

7 

9,541 

48 

6 

3 

7.358 

98 

13 

3 

9,541 

49 
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2 

7.358 
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1 
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TABLE  4 < CONT » ) 


ALPHA=  .05  Al=l/3  K=2 
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1 

0 

0 
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1 

0 

2.125 
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3.162 
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10 
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4.341 
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16 
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5.168 
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5.168 

20 
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21 
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24 
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26 
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28 

5 
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10 

0 

6.102 

30 

3 

3 
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31 

8 

1 
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32 

1 

4 

6.241 

33 

6 

2 

6.241 

34 

11 

0 

6.241 

35 

4 

3 

6.288 

36 

9 

1 

6.359 

37 

2 

4 

6.548 

38 

0 

5 

7.009 

39 

7 

2 

7.009 

40 

12 

0 

7.009 

41 

5 

3 

7.009 

42 

10 

1 

7,009 

43 

3 

4 

7.009 

44 

1 

5 

7.148 

45 

8 

2 

7.148 

46 

13 

0 

7,148 

47 

6 

3 

7.148 

48 

4 

4 

7.274 

49 

11 
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2 
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7.457 
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51 

9 

2 
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52 

0 
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7.894 
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14 

0 

7.894 
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7 

3 

7.894 

55 

5 

4 

7.894 
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12 
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7.894 

57 

3 

5 

7.894 

58 

10 

2 

7.894 

59 

1 

6 

8.035 

60 

8 

3 

8.035 

61 

15 

0 

8.035 

62 

6 

4 

8.036 

63 

13 

1 

8.036 

64 

4 

5 

8.175 

65 

11 

2 

8.177 

66 

2 

6 

8.344 

67 

9 

3 

8.345 

68 

16 

0 

8.345 

69 

0 

7 

8.765 

70 

7 

4 

8.765 

71 

14 

1 

8.765 

72 

5 

5 

8.765 

73 

12 

2 

8.765 

74 

3 

6 

8.765 

75 

10 

3 

8.765 

76 

17 

0 

8.765 

77 

1 

7 

8.906 

78 

8 

4 

8.906 

79 

15 

1 
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80 

6 

5 

8.907 

81 

13 

*? 

8.915 

82 

4 

6 

9.059 

S3 

11 

3 

9.061 

84 

18 

0 

9.061 

85 

o 

Am 

7 

9.134 

86 

9 

4 

9.134 

87 

0 

8 

9.142 

88 

16 

1 

9.144 

89 

7 

5 

9.249 

90 

14 

Atm 

9,287 

91 

5 

6 

9.344 

92 

12 

3 

9.344 

93 

3 

7 

9.568 

94 

19 

0 

9.568 

95 

10 

4 

9.631 

96 

1 

8 

9.631 

97 

17 

1 

9.631 

98 

8 

5 

9.631 

99 

15 

2 

9.765 

100 

6 

6 

9.765 
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TABLE  4 ( CQNT » ) 


ALPHA=  « 05  Al= . 40  K=2 
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THETA* 

1 

0 

0 

1.797 

2 

1 

0 

2.056 

3 

0 

1 
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4 

2 

0 
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5 

1 

1 
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6 

3 

0 

3.158 

7 

0 

o 

Am 

3.777 

8 

2 

1 
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9 

4 

0 

3.777 

10 
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4.004 
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0 
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0 
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Am 
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15 

4 

1 
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16 

6 
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1 
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4.876 
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3 

2 

4.985 

19 

5 

1 

5.086 

20 

0 

4 

5.091 

21 

7 

0 

5.206 

22 
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«• 

3 

5.206 

23 

4 

2 

5.464 

24 

1 

4 

5.713 

25 

6 

1 

5.713 

26 

3 

3 

5.831 

27 

8 

0 

5,831 

28 

0 

5 

5.933 

29 

5 

2 

5.933 

30 

2 

4 

6,108 

31 

7 

1 

6,130 

32 

4 

3 

6.213 

33 

9 

0 

6.213 

34 

1 

5 

6.527 

35 

6 

2 

6.528 

36 

3 

4 

6.650 

37 

8 

1 

6,651 

38 

0 

6 

6.774 

39 

5 

3 
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40 

10 

0 

6.775 

41 

'y 

A. 

5 

6.903 

42 

7 

o 

An 

6.914 

43 

4 

4 

7.113 

44 

9 

1 

7.151 

45 

1 

6 

7.324 

46 

6 

3 

7.324 

47 

11 

0 

7.324 

48 

3 

5 

7.451 

49 

8 

2 
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50 

0 

7 
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I 

XI 
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51 

5 

4 
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52 

10 

1 
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53 

2 

6 
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54 

7 

3 

7.724 

55 

12 

0 

7.724 

56 

4 

5 

7.908 

57 

9 

o 

Am 

7.911 

58 

1 

7 

8.106 

59 

6 

4 

8.107 

60 

11 

1 

8.107 

61 

3 

6 

8.237 

62 

8 

3 

8.238 

63 

13 

0 

8.238 

64 

0 

8 

8.259 

65 

5 

5 

8.371 

66 

10 

2 

8.375 

67 

2 

7 

8.572 

68 

7 

4 

8.573 

69 

12 

1 

8.573 

70 

4 

6 

8.579 

71 

9 

3 

8.579 

72 

14 

0 

8.579 

73 

1 

8 

8.878 

74 

6 

5 

8.878 

75 

11 

2 

8.878 

76 

3 

7 

9.010 

77 

8 

4 

9.013 

78 

0 

9 

9.065 

79 

13 

1 

9.074 

80 

5 

6 

9.114 

81 

10 

3 

9,155 

82 

2 

8 

9.339 

83 

15 

0 

9.339 

84 

7 

5 

9,339 

85 

12 

o 

Am 

9.339 

86 

4 

7 

9.344 

87 

9 

4 

9.457 

88 

1 

9 

9,461 

89 

14 

1 

9.461 

90 

6 

6 

9,591 

91 

11 

3 

9,640 

92 

3 

8 

9.775 

93 

16 

0 

9.775 

94 

8 

5 

9.778 

95 

0 

10 

9.814 

96 

13 

2 
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97 

5 

7 
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98 

10 

4 
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99 

2 

9 

10.095 
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15 

1 
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